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1 Introduction 



Among the most attractive options of facilities at the high energy frontier of elementary 
particle physics are high luminosity e'^e~ linear colliders like TESLA the NLC 

or the JLC f^. However, going to higher energies and higher luminosity becomes a real 
challenge for working out Standard Model (SM) predictions of the adequate precision because 
of the dramatically increasing complexity of perturbative calculations. Here we consider all 
four fermion production processes in electron positron annihilation together with real photon 
emission at the tree level. These processes with 6 and 7 external particles at the tree level 
are described by from 10 up to 1008 Feynman diagrams in a given channel, neglecting the 
Higgs boson coupling to light fermion flavors, and the physical cross section is the result 
of a very obscure quantum mechanical interference between all these diagrams. Interesting 
are of course those cases where the result is dominated by a few diagrams like in VT-pair 
production where we have three relevant "signal diagrams". However, also in these cases 
which allow a relatively simple physical interpretation, many other diagrams may play a role 
as a background contribution which affects the precise interpretation of the signal process. 
The latter give us information about the gauge boson parameters and the triple or quadruple 
gauge couplings. Similarly, the properties of the Higgs boson will be fixed by its contribution 
as unstable intermediate state. The most important physics cases have been reviewed in [Q, 
for example. A number of dominant e~^e~ 4f channels have been explored experimentally 
at LEP-II (1996-2000) |^ and the measurements confirmed SM predictions at the level of the 
most relevant 0{a) corrections. 

In the approximation of massless fermions all possible four fermion channels e'^e" —>■ 4f have 
been investigated in (EXCALIBUR) and those for e~^e~ —>■ Af'j in (RacoonWW). Here we 
extend these investigations to a calculation with nonzero fermion masses. Keeping nonzero 
fermion masses will be important in cases where predictions at the 1% accuracy level are 
required [§, ^. Finite masses also provide a natural regularization of distributions which 
become singular in the massless limit. Massive calculations thus provide reliable benchmarks 
for massless calculations with cuts. The latter are much simpler and hence much faster than 
calculations with massive codes. The hard bremsstrahlung processes are of interest in their own 
right and may be used to investigate anomalous WW'j and WW'yy couplings, for example. 

Our calculation is considered to be a building block (the soft plus hard bremsstrahlung part) 
for a complete 0{a) calculation of the processes e~^e~ —>■ 4f. Such calculations have been 
attempted in [|1^ (see also lll|). This would also extend existing calculations of I^-pair 



production in the double pole approximation [|l^] (RacoonWW) and [0 (KORALW/YFSWW) (see 



also [|T^ (EEWW)) which incorporate the one- loop corrections for production of on-shell Im- 



pairs JTSf and their subsequent decay into fermion-pairs [jT6|. 



There already exist a number of codes which allow to calculate exact matrix elements for 
e~^e~ — i> 4f, 4f7 for massive fermions. Some of the program packages available are general 
purpose packages which allow for an automatic calculation of tree-level amplitudes and for 
their numerical evaluation. Known programs, which may be utilized for tree-level calculations 



of the kind we are interested in are: GRACE/BASES MADGRAPH/HELAS ]T8[, CompHEP [jT9 
(squared matrix element technique), WPHACT |20|, NEXTCALIBUR ||2ll (initial state radiation 
photons generated via the structure function approach), HELAC/PHEGAS (recursive Dyson- 



Schwinger equation approach), WRAP |2§ (ALPHA dJ algorithm) and 'Mega/WHIZARD p| 
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Most of the codes work on the basis of hehcity amphtudes and some use the structure function 
approach to generate the photons. Except for NEXTCALIBUR, which is speciahzed to e~^e~ — > 4f , 
all other programs allow to generate and evaluate amplitudes for other type of processes. For 
more details and comparisons we refer to the Four Fermion Working Group Report ||2^ . 



Many channels have their own specific problems concerning numerical stability and/or effi- 
ciency and need separate consideration and optimization. We therefore present, in this paper, 
a different approach which is optimized for each individual channel. We present reference ta- 
bles for cross section calculations at y/s = 200 GeV and 500 GeV for all leptonic, semi-leptonic 
and hadronic channels. The present work completes previous investigations of specific channels 
presented in P, 0. 

Our calculation has to be extended to include the 0{a) virtual corrections in future. Precise 
knowledge of the various channels is crucial for the precise determination of properties of the 
unstable gauge and Higgs bosons as well as to reveal possible anomalous couphng ^ |28 
which might exist beyond the SM. 

In the following we outline our calculation, present the numerical results and end with the 
conclusions. 



2 Calculation 



The matrix elements of the reactions 

e+e- ^ 4/ (1) 

and 

e+e- ^ 4/7 (2) 

are calculated by utilizing the helicity amplitude method described in p. As in g], the 
photon propagator is taken in the Feynman gauge while for the propagators of the massive 
gauge bosons we use the unitary gauge. Constant widths of the electroweak gauge bosons, 
T]y,Tz, Higgs boson, Th, and the top quark, are introduced through the complex mass 
parameters 

My = my - imvTv, V = W,Z, = m]j - iniHT h, Mt = rrit - iTt/2, (3) 

in the corresponding propagators 



q^-M^y ' ^^^^^ - 'J^' ^^^^^ - 



both in the s- and t-channel Feynman diagrams. The electroweak mixing parameter is kept 
real 

sin^ 6*1^ = 1 — "m^/ml. (5) 

This kind of parametrization is usually referred to as the fixed- width scheme (FWS). Our 
results presented in the next section have been obtained in the FWS. In our program ee4f7, 
it is also possible to define sin^ 6^ in terms of the complex masses of (^) as 

sin2^H^ = l-M^/M|, (6) 
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which is usually called the complex- mass scheme (CMS). The CMS has the advantage that 
it satisfies the SU{2) x f/(l) Ward identities at tree level 0. However, the fact that (|) 
makes some of the SM couplings complex quantities may become a source of discomfort. In 
the FWS on the other hand, all the couplings remain real. Only the electromagnetic gauge 
invariance is satisfied exactly, however, and this only provided that Tt and the other fermion 
widths are vanishing. It should be stressed at this point, that for vanishing fermion widths, 
electromagnetic gauge invariance is preserved with non-zero fermion masses and with the 
gauge boson widths Tw and Tz treated as independent parameters. If a non-vanishing top 
quark width is introduced through substitution (|^), which is done in order to regularize the 
on-mass-shell pole of a top quark propagator at tree level in reactions (|^) containing a single 
top quark in the final state, the external electromagnetic gauge invariance gets violated. It 
can be restored by redefining the Dirac bi-spinor representing the external top quark in such 
a way that it satisfies the Dirac equation with the complex top quark mass of Eqs. @, which 
would obviously require a complex top quark four momentum in the phase space generation. 
As the production of a on-mass-shell top quark is a rather unphysical process, one should not 
be to much concerned about the problem. Fortunately, as we shall see in the next section, the 
violation of the gauge symmetry does not lead to dramatic effects for the total cross sections 
of reactions (§) containing a top quark in the final state. A more realistic treatment of the 
top quark has to include its decay and thus requires the consideration of e~^e~ — >• 6/, 6/7 
channels, which is beyond the task of the present investigation. 

The hadronic channels are discussed only at the level of quark parton production. Quark- 
mass effects will be estimated by adopting the so-called current-quark masses in the MS 
scheme at a scale /i ~ 2 GeV. This should allow us to get an idea about the size of mass 
effects and eventually allow us to establish suitable cuts which eliminate the mass sensitivity 
of quark production cross sections. In any case, taking into account mass effects, provides an 
improvement over calculations in the approximation of massless quarks. For observables which 
exhibit a substantial mass dependence of course one would have to discuss more carefully the 
precise physical meaning of quark masses in the given process. 



3 Results 



In this section, we will present numerical results for all the four-fermion channels of reactions 
(|I|) and (0) which are possible in the SM. 

We define the SM physical parameters in terms of the gauge boson masses and widths, the top 
mass and width, and the Fermi coupling constant. We take the actual values of the parameters 
from [p9|] : 

= 80.419 GeV, Fvf = 2.12 GeV, = 91.1882 GeV, F^ = 2.4952 GeV, 

mf = 174.3 GeV, = 1.16639 x 10"^ GeV'^. (7) 

We assume a Higgs boson mass of mj^^ = 115 GeV and calculate the Higgs boson width with 
the lowest order formula in the SM. The top quark width is assumed to be F^ = 1.5 GeV. 
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For the sake of definiteness we also list the other fermion masses that we use in the calcula- 
tion pi: 



me = 0.510998902 MeV, = 105.658357 MeV, = 1777.03 MeV, 

mu = 5 MeV, nia = 9 MeV, m, = 150 MeV, = 1.3 GeV, rrib = 4.4 GeV. (8) 

We neglect the Cabibbo-Kobayashi-Maskawa (CKM) mixing, i.e., we assume the CKM matrix 
to be the unit matrix. However, it is possible to run the program with nontrivial CKM mixing 
as well. 

The effective fine structure constant (at scale ~ M\y) is calculated via 

aw = V2Gfj^mly sin^ 6w/ tt (9) 

utilizing the real electroweak mixing parameter sin^^vF defined by (^. In ee4f7, it is also 
possible to perform computations with the complex sin^ 6iy of Eq. (Mj and the complex 



w 



of Eqs. (g), i.e. with the complex aw- The photon coupling to fermions and gauge bosons 
is given by the fine structure constant in the Thomson limit a = 1/137.0359895 and the 
quark-gluon strong interaction "constant" by as{Mz) = 0.1185. 



We will apply the following set of "standard cuts" which have been proposed in [2q: 



cos 9{l, beam) < 0.985, ^(7, /) > 5°, > 1 GeV, m(g, q') > 10 GeV , 

cos ^(7, beam) < 0.985, ^(7,9) > 5°, > 5 GeV, (10) 

where I, q, 7, and "beam" denote charged leptons, quarks, photons, and the beam (electrons 
or positrons), respectively, and 0{i,j) the angles between the particles i and j in the center of 
mass system. Furthermore, m{q, q') denotes the invariant mass of a quark pair qq' . Note that 
we are not applying a corresponding cut to the invariant mass of the charged lepton pairs. 

The errors we will quote in the Tables below have been evaluated as follows: For each separate 
channel of the multichannel Monte Carlo (MC) integration the error is calculated by VEGAS [|| 



This is a purely statistical error equivalent to one standard deviation. We added linearly 
standard deviations for all the channels used in an integration and this is what is our error. 
This provides a more conservative estimate for the error than for example adding up partial 
errors in quadrature. 

Except for the check of electromagnetic gauge invariance, discussed in the previous section, 
we perform a few other checks. Whenever the fermion masses play no role we have reproduced 
the results of 0. The matrix elements of almost all channels of the processes (|l|) and (0) under 



consideration have been checked against MADGRAPH [|18|. The comparison was not simple for 
the channels involving a gluon exchange, since the version of MADGRAPH which we are using, 
generates either the electroweak or the QCD part, but not both simultaneously. In addition, 
for the reaction e^e~ — > e'^e~e"'"e~7 MADGRAPH generates only 999 instead of all 1008 Feynman 
graphs. The phase space generation routines have been thoroughly checked against each other 
before they have been combined into a multichannel phase space generation routine. The total 
cross sections of the reactions (P containing a single top quark in the final state 

e+e- tiff', (11) 
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Table 1: Cross sections in femto-barns (1 fb=10~"'^^ barns) of e^e~ — > csr~i/^7 and e'^e" thT~i>T-^ 
in CMS at ^/s = 500 GeV for different photon energy cuts lo. as, corresponding to -E^ < uj, is the 
cross section in the soft photon hmit and a^, corresponding to E.y > uj, is the hard bremsstrahlung 
cross section. A fictitious photon mass m-y = 10~^ GeV has been introduced in order to regularize 
the infrared divergence. No other cuts except for uj and are present. 



UJ 

(GeV) 






7 / ^/ 








as (fb) 


OK (fb) 


Os + Oh (fb) 


Os (fb) 


Oh (fb) 


Us + Oh (fb) 


0.1 


186.37(9) 


250.7(3) 


437.1 


51.85(3) 


56.88(9) 


108.7 


0.01 


108.89(6) 


328.1(3) 


437.0 


32.29(2) 


76.1(1) 


108.4 


0.001 


31.39(2) 


405.5(4) 


436.9 


12.729(8) 


95.3(1) 


108.0 



where / = e , /i , r , li, s and /' = z/g, z/^, z/^, m, c, respectively, have been calculated in an ar- 



bitrary linear gauge This does not allow to estimate the absolute size of gauge violation 
effects caused by the nonzero widths of the unstable fermions. However, as the transition be- 
tween two linear gauges, the 't Hooft-Feynman and unitary gauge, which has been numerically 
performed by changing the gauge parameter from 1 to 10^^ has caused practically negligible 
change in the cross sections at center of mass energies up to 2 TeV, typical for a linear collider, 
one may expect that the gauge violation effects are not very dramatic for total cross sections. 

Moreover, another test of reliability of our results has been performed for several final states. 
We have split the cross section of the bremsstrahlung process (0) into a soft photon part 
(Ts, which includes the photons with energies < u, and a hard photon part ah, including 
the contributions from photon with energies > u, and checked whether the combined 
bremsstrahlung cross section = Us + Oh is independent of the photon energy cut p. In 
Table 1, we illustrate this independence, which holds within one standard deviation of the MC 
integration, for e+e^ —>■ cst~Ut-'j in the CMS. At the same time, in Table 1, we see a small 
dependence on the cut-off parameter uj, which is at the level of about two standard deviations, 
for the bremsstrahlung reaction e~^e~ —>■ thr'Vrl- The cut dependence is most probably 
induced by the violation of external electromagnetic gauge invariance caused by the nonzero 
top quark width, which has been introduced to e^e~ tbr'Ur'y in a somewhat asymmetric 
way, related to the fact that the top quark is regarded as on-mass-shell particle at the same 
time when the anti-top quark decays. It should be stressed, that the soft bremsstrahlung 
cross sections as presented in Table 1 are unphysical, as they contain the unphysical photon 
mass m^. They are only given in order to show that the leading logarithmic contributions are 
treated properly within the CMS for reactions which do not contain nonzero fermion width 
and to illustrate the size of cut off dependence caused by the nonzero top quark width. As the 
cut dependence in e~^e~ — >■ tbr'h'r'y is of the order of 1% of the corresponding four fermion 
Born cross section, one should certainly elaborate more on this issue in future, in the context 
of a more realistic six fermion reactions, which would treat the decay of the top quark on the 
same footing as that of the anti-top quark. 

We compare our results for the total cross sections of e+e^ cs/i~z/^ at ^/s = 200 GeV with 
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Table 2: Cross sections in fb of e^e cs/j, at \/s = 200 GeV for different cuts on the photon 



angle with respect to the quarks 6{'j,q) or muon 6(^,11) and the remaining cuts as in (10). Here we 
use physical parameters of |23] and parametrize the photon couplings by a = 1/137.0359895. 







[23| 


Present work 


5° 


5° 


74.294(29) 


74.267(60) 


1° 


1° 


93.764(37) 


93.70(7) 


5° 


r 


90.157(36) 


90.13(7) 


5° 


0.1° 


104.777(46) 


104.78(7) 


5° 


0° 


105.438(45) 


105.48(7) 



Table 3: e+e — > 4f cross sections a and e+e Afj cross sections in fb at t/s = 200 GeV and 
^/s = 500 GeV for different four-fermion final states corresponding to the Ty"'"M^~-pair signal. The 



cuts are those of (IC) 



Final 

state 


v/i = 200 GeV 


= 500 GeV 







a 






630.65(31) 


70.547(83) 


211.11(13) 


23.601(46) 


udT'Vr 


630.18(31) 


68.321(74) 


210.95(13) 


23.386(44) 




630.40(31) 


69.501(80) 


211.03(13) 


23.285(47) 




629.93(31) 


67.279(72) 


210.87(13) 


23.077(43) 








58.88(30) 


9.467(60) 








58.80(29) 


9.295(56) 


csdu 


1838.6(1.4) 


172.74(28) 


749.07(50) 


68.34(23) 


tbdu 






177.8(1.9) 


25.55(42) 


tbsc 






177.4(1.9) 


25.37(37) 




205.88(15) 


25.784(44) 


60.762(62) 


7.842(23) 


uddu 


1921.4(7) 


188.19(46) 


780.66(25) 


74.99(28) 


esse 


1925.7(8) 


184.07(46) 


782.62(28) 


73.46(25) 


tbbi 






0.85519(56) 


0.073748(78) 




218.91(19) 


28.232(55) 


63.933(70) 


8.475(26) 




214.94(20) 


26.280(50) 


63.468(74) 


8.299(20) 


Vf,e^e~V(, 


259.55(31) 


32.012(93) 


195.22(42) 


24.85(14) 
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Table 4: e~^e — > 4f, 4f7 cross sections in fb at ^/s = 200 GeV and ^/s = 500 GeV for different 
four-fermion final states corresponding to the single- signal. The cuts are those of ([lO|). 



Final 

state 


v/i = 200 GeV 


^ = 500 GeV 


cr 


dry 


a 


CT-y 


ude~Ue 


661.68(40) 


72.95(10) 


354.06(27) 


38.876(88) 




661.42(40) 


71.84(10) 


353.91(27) 


38.371(88) 








58.54(65) 


9.43(16) 




216.29(21) 


27.473(57) 


107.34(14) 


13.677(43) 




216.13(21) 


26.709(55) 


107.25(13) 


13.471(42) 



those of 1^ in Table 2. As in different cuts on the photon angle with respect to the 
quarks 0{'-f,q) or muon 6{'~f,fi) are imposed while the remaining cuts are those given in (0). 
For the comparison we use the physical parameters of i.e. = 1.16637 x 10~^ GeV~^, 



w 



mz = 91.1867 GeV, mw = 80.35 GeV, sin^ Ow = 1 - Tz = 2.49471 GeV, F 

2.04277 GeV, = 0.10565839 GeV, m, = 0.15 GeV, = 1.55 GeV. Although it has 
not been explicitly specified, we assume that is using a = 1/137.0359895 for the photon 



coupling strength. We find that the results agree perfectly within one standard deviation of 
the MC integration. 

Our results for all channels of reactions e~^e~ 4f, 4f7 possible in the SM are collected 
in Tables 3-7. We present total cross sections at two center of mass energies, = 200 
GeV and a/s = 500 GeV, with cuts defined by (|TD|), except for e~^e~ e~^e~e~^e~, where 
we have imposed another cut on the angle between the final state electrons and/or positrons, 
^(e±,e±) > 5°. 

In Tables 3 and 4, we show the results for the channels corresponding to the iy"'"iy~-pair 
and single-iy signal. These channels are usually classified as charged current reactions. The 
relative magnitude of the cross sections in both tables reflects the naive counting of the color 
degrees of freedom, e.g., the cross sections of purely hadronic channels are about a factor 3 
bigger than the cross sections of semi-leptonic channels and the latter are a factor 3 bigger 
than the cross sections of purely leptonic reactions. This somewhat general rule is obviously 
violated in reactions which receive contributions from the gluon or t-channel photon and Z 
exchange. Except for the channels containing heavy quarks, t and b, for lighter flavors, the 
fermion mass effects are not big. However, for individual channels they are of the order of 
a few per cent, as it has been already pointed out in It is amazing that the mass effect 
is inverse for e~^e~ uddu and e^e~ esse. The inversion is not due to the Higgs boson 
exchange, but in fact is a consequence of the cuts ( p!0D which we imposed. The latter reduce 
the contribution of the s-channel Feynman diagrams to e+e^ uddu to much larger extent 
than to e^e~ — esse, because the cuts on the invariant mass of the quark pairs restrict the 
phase space much more severely for lighter quarks than for heavier ones. Without the cuts, 
the cross section of e"'"e~ uddu becomes bigger than that of e'*'e^ CSSCj cLS expected. We 
do not show cross sections at ^/s = 200 GeV for reactions containing a t-quark in the final 
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states, as they are negligibly small 



m 



Table 5: Cross sections in fb of the purely leptonic neutral-current channels of (|T|) and (|2|) at ^/s = 200 
GeV and ^/s = 500 GeV. The cuts are given by ( p^ ) with the exception of e'^e^ — > e'^e~e~^e~, 
where we have imposed another cut on the angle between the final state electrons and/or positrons, 
e(e±,e±) > 5°. 



Final 

state 


= 200 GeV 


^ = 500 GeV 


a 


/ 


a 


7 




10.267(14) 


2.1787(91) 


2.5117(44) 


0.6495(40) 




12.729(10) 


1.6998(44) 


3.0336(31) 


0.5695(20) 




9.1659(59) 


1.2307(27) 


2.7174(27) 


0.5221(15) 




10.608(6) 


0.57713(82) 


3.9179(43) 


0.4645(11) 




137.18(90) 


12.93(31) 


43.80(38) 


4.58(12) 




54.49(18) 


7.115(55) 


16.860(61) 


2.685(23) 




17.780(18) 


2.1467(43) 


24.031(92) 


3.479(20) 




11.721(10) 


1.4797(24) 


21.389(90) 


3.432(21) 




11.448(8) 


0.6033(7) 


24.728(47) 


2.349(7) 




23.503(22) 


2.7789(68) 


9.453(14) 


1.3076(57) 




6.747(13) 


1.4307(84) 


1.5106(36) 


0.3860(30) 




3.7283(30) 


0.7943(21) 


1.0341(11) 


0.2727(8) 




5.2660(23) 


0.28562(25) 


1.9577(15) 


0.23149(34) 




50.53(10) 


5.770(58) 


13.927(32) 


2.163(21) 




5.9815(27) 


0.30563(24) 


22.482(61) 


2.0856(92) 



The results for the neutral current channels of reactions (|1]) and (^ are shown in Tables 5, 
6 and 7. We list purely leptonic channels in Table 5, semi- leptonic channels in Table 6 and 
purely hadronic channels in Table 7. The cross sections in Tables 5, 6 and 7 are typically much 
smaller then those of Tables 3 and 4. Mass effects on the other hand are bigger. The stronger 
dependence on fermion masses can be explained as follows. The neutral-current reactions are 
dominated by s-channel Feynman diagrams which contain the propagator of a photon decaying 
into a fermion pair. This causes a ~ V'^//' behavior of the matrix element squared and results 
in a relatively high sensitivity to the fermion pair threshold Sff = {rrif + mf)'^. There is a 
relatively big effect for charged lepton pairs t~^t~ and much smaller effect for quark 

pairs, except for tt of course. This is due to the fact that there is no cut on the invariant 
mass of a charged lepton pair in Eqs. ([101) . Again we observe an inverse mass effect in some 
channels, especially those containing a neutrino pair. The cross section of e+e^ —>■ ddVf.Vf, is 
bigger than that of e^e~ — > uuVeVe although the mass of the cZ-quark is almost twice as big as 
that of the w-quark. The inverse mass effect is caused by the invariant mass cut m(g, g') > 10 
GeV of ([10|), which is more restrictive for lighter fermion pairs than for heavier ones and by 
the fact that there is neither an invariant mass cut nor an angular cut on a neutrino-pair. We 
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Table 6: Cross sections in fb of the semi-leptonic neutral-current channels of ([T|) and (|2|) at ^/s = 200 
GeV and = 500 GeV. The cuts are those specified by (|lO[). 



Final 


= 200 GeV 


^ = 500 GeV 


state 


cr 




cr 






27.341(30) 


4.874(13) 


6.855(10) 


1.5630(62) 




19.543(17) 


3.4694(98) 


5.9096(69) 


1.3614(49) 




22.614(14) 


2.2049(34) 


8.3860(99) 


1.2897(30) 




27.287(27) 


4.826(19) 


6.9946(88) 


1.5924(74) 




19.560(16) 


3.4407(92) 


6.0566(68) 


1.3922(48) 




22.655(13) 


2.1695(44) 


8.6652(83) 


1.3434(39) 








0.1832(2) 


0.02165(5) 


ttr+T~ 






0.11991(10) 


0.01554(2) 








0.07010(2) 


0.004577(4) 




29.541(24) 


4.229(13) 


7.0683(80) 


1.3715(60) 


ddr^T^ 


21.256(16) 


3.0660(51) 


6.2854(61) 


1.2494(26) 


ddv^v^ 


24.634(15) 


1.5960(28) 


9.0094(91) 


1.1230(33) 




29.542(25) 


4.240(14) 


7.0702(73) 


1.3737(55) 




21.257(16) 


3.0666(76) 


6.2907(62) 


1.2489(38) 




24.637(15) 


1.5975(29) 


9.0125(89) 


1.1301(29) 




29.536(23) 


4.150(13) 


8.6899(74) 


1.7290(56) 


bbr^T^ 


21.349(16) 


3.0179(76) 


7.9184(66) 


1.5972(43) 


bbu^u^ 


25.019(15) 


1.5457(28) 


12.322(9) 


1.5911(34) 




77.31(26) 


10.438(74) 


28.37(11) 


4.432(33) 




77.02(19) 


10.321(61) 


29.030(91) 


4.532(32) 


ite+e- 






172.63(30) 


15.657(45) 


dde^e^ 


58.759(82) 


7.308(27) 


23.668(46) 


3.363(21) 




58.734(80) 


7.337(26) 


23.702(46) 


3.377(18) 


bbe^e^ 


58.082(66) 


7.125(21) 


29.38(11) 


4.152(27) 




25.502(19) 


2.4289(41) 


49.18(22) 


6.300(38) 




25.907(21) 


2.4210(38) 


55.13(26) 


7.180(47) 








0.07400(4) 


0.004810(4) 


ddugUg 


26.820(20) 


1.6916(23) 


57.03(27) 


5.731(36) 




26.824(22) 


1.6955(24) 


57.320(28) 


5.777(41) 


bbUel^e 


31.346(28) 


1.8819(33) 


128.31(65) 


12.59(12) 
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Table 7: Cross sections in fb of the purely hadronic neutral-current channels of (||) and (|2|) at 
^ = 200 GeV and ^/s = 500 GeV. The cuts are those specified by (0). 



Final 
state 


v/i = 200 GeV 


^ = 500 GeV 


cr 


cr^ 


cr 


cr^ 


uuss 


83.63(12) 


13.45(13) 


32.561(51) 


5.891(41) 


uubb 


86.565(95) 


14.39(11) 


38.475(45) 


7.004(36) 


ccdd 


84.41(12) 


13.742(76) 


33.150(92) 


5.969(41) 


ccbb 


87.209(83) 


14.445(99) 


39.237(42) 


7.149(35) 


itdd 


— 


— 


0.903(1) 


0.0801(2) 


itss 


— 


— 


0.903(1) 


0.0803(2) 


ddss 


79.658(95) 


9.77(11) 


29.812(39) 


4.159(34) 


ddbb 


82.486(87) 


10.436(86) 


37.435(39) 


5.363(26) 


uucc 


87.06(13) 


15.98(12) 


35.826(61) 


7.581(52) 


uutt 


— 


— 


0.9071(11) 


0.09440(21) 


uuuu 


42.465(45) 


7.794(35) 


17.410(27) 


3.684(19) 


cccc 


43.323(44) 


7.935(27) 


18.076(23) 


3.867(14) 


dddd 


39.173(32) 


4.791(20) 


14.758(15) 


2.0706(76) 


ssss 


39.167(32) 


4.810(20) 


14.767(15) 


2.0769(76) 


bbbb 


41.667(24) 


5.217(14) 


22.242(14) 


3.2142(67) 



observe that the mass effects depend on cuts. They may be quite different for different choices 
of cuts. 

Whether or not the mass effects will play a role in the analysis of the future data depends 
mostly on the luminosity of a future linear collider. If we assume a total integrated luminosity 
of 500 fb~^ they will certainly become relevant. Therefore it is better to keep nonzero fermion 
masses in the calculation, or test the massless fermion generators against a massive one for 
each given set of kinematical cuts. 

If we compare the channels containing a bb- or cc-pair with the channels containing lighter 
quark flavors, we see a clear signal of the Higgs-strahlung reaction e"'"e~ ZH, especially at 
a/s = 500 GeV, which exceeds the ZH production threshold for a Higgs boson mass of 115 
GeV. In case of e"'"e~ bbuei^e, ccue^e, we see also a signal of a W^W~ fusion mechanism 
of the Higgs boson production. For the final state containing a 66-pair, the signal is visible 
already at ^/s = 200 GeV and it becomes much more pronounced at ^/s = 500 GeV. 

As the threshold energy for e^e^ — > itit is bigger than 500 GeV, we do not show its cross 
section in Table 7. However, even at ^/s = 800 GeV the cross section of e+e~ — > itit without 
cuts is of the order of 10"^ fb. 
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4 Conclusions 



We have developed a program package ee4f 7 which allows us to calculate in an efficient manner 
any process e"'"e~ — >■ 4f , 4f7 keeping nonzero fermion masses. Corresponding results have been 
investigated and discussed for all channels at reference energies ^/s = 200 GeV and 500 GeV. 
For detailed investigations of physics at a high luminosity linear collider like TESLA, these 
mass effects should be taken into account. Nonzero fermion masses also provide a physical 
regularization of matrix elements which exhibit collinear singularities in the zero mass limit. 
The results thus may be used as benchmarks for calculations performed in the approximation 
of vanishing fermion masses which allows one to perform calculations with much less numerical 
efforts. Our results may be considered as part of the complete 0{a) corrections which we will 
need to understand physics at future high energy linear colliders in a precise manner. The 
calculation of the missing virtual corrections to e~^e~ 4f is one of the big challenges for the 
future. 
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